General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



(NASA-OR-152187) I MPLENENTATION OP AN 
OPTIMUM PROFILE GUIDANCE SYSTEM ON STOLAND 
(Analytical Mechanics Associates, Inc.) 

49 p HC A03/MP A01 CSCL 17G 

G3/04 


N79-10038 


Tnclas 

36091 



ANALYTICAL MECHANICS ASSOCIATES, INC. 


24H \ OLD Ml DDL! I IHIHUV 
MOON I MN VILW. CALIFORNIA 94041 


NASA CH-152187 


IMPLEMENTATION OL AN OPTIMUM 
PROFILE GUIDANCE SYSTEM ON 

STOLAND 


By 

Paul F. Flanagan 


September 19?* 


Distribution of this, report is provided in the interest of information 
exchange. Responsibility for the contents resides 
in the author or organization that prepared it. 


Prepared under Contract No. NAS2-9460 
By 

Analytical Mechanics Associates, Inc. 

Mountain View , California 

for 

AMES RESEARCh CENTER 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Moffett Field, California 



ABS1RAC! SIMMARY 


Thi report describes briefl) the implementation on the STOLAND airborne 
digital computer of an Optimum Profile guidance system for the Augmentor Wing 
Jet STOL Research Aircraft. Major tasks under this contract were to implement 
the guidance and control logic, developed by NASA, to airborne computer software 
and to integrate the modul' with the existing STOLAND navigation, display, and 
autopilot routines. The optimum profile guidance system comprises an algorithm 
for synthesizing minimum fuel trajectories for a wide range of starting positions 
in the terminal area and a control law lor flying the aircraft automatically along 
the trajectory. The report also touches on operational aspects of the system. 

In addition to describing the avionics software developed, a FORTRAN program 
is described that has been constructed to reflect the modular nature and algorithms 
implemented in the avionics software. The technical monitors as well as the 
principal contributors to the analytical development of the system are John D. 
McLean and Heinz Erzbergcr of NASA/ARC. 
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1 PROJEC1 SYSTEM DESCRIPTION 


Background 

The NASA Ames STOLANl) system is an Integrated digital avionics package 
designed for testing guidance, navigation and control concepts and for investigating 
operational procedures for short takeoff and landing aircraft. The STOLAND 
system includes njvaid receivers, onboard sensors and pilot control and command 
inputs that are interfaced with a Sperry 1&19A computer. The computer is used 
to provide both pilot assist modes and automatic modes for control of the aircraft. 
The Sioland system provides various displays (EADI, MFD, I1S1 > , control actuators, 
mode select and data entry devices, vehicle sensors and a data acquisition system. 
The Sperry 1819A computer with auxiliary memory used in the STOLAND system 
provides 32K of words of lto-bit memory’. 


The standard Sperry developed software for the At gmentor Wind aircraft 
(Reference 1) provides interface with all requit vd systems to support flight opera- 
tions. In the implementation of the Optimum profile guidance ana control program 

significant revision of the standard software svstem was required to pro- 
vide required memory and CPI’ time. The following revisions to the standard 

system were performed. 

a) removal of the existing reference flight path package to create memory 
space 

b) implement a revised executive package to provide sufficient real-time 
for execution of the Optimum profile guidance and control pack;tge 

c) develop a revised utility package with overlay capability to create memory 
space and provide a more comprehensive diagnosis tool for laboratory 
checkout (Reference 2) 

d) Revise the strip chart routine to provide multiplexed strip charts in the 
laboratory to provide an improved diagnosis tool 
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e) moving all code (rom bank 5 to provide an effective data area for the 
Optimum profile guidance and control package 

f) revision of the Keyboard entr> and display MFD, KADI, anti HSI routines 
lo provide the required interface 

Design Method for the Guidance System 

The purpose of the system is to assess the performance of a fuel efficient 
onboard guidance system for powered lift of STOL aircraft. The system generates a 
minimum fuel noise trajectory for the Augmcntor Wing aircraft and controls the 
aircraft along that trajectory. 

Tpon entering the terminal area the system, operating in fast-time, 
synthesizes a minimum-fuel reference trajectory from the current aircraft position 
to touchdown. The calculation makes ise of energy -rate schedules, generated off 
line and stored in the STOLAND computer. The schedules give the controls (flaps, 
throttle, noz/lc, and angle of attack) yielding the least fuel flow for a given normal- 
ized energy rate, E , aircraft weight, altitude and equivalent airspeed. The 

normalized energy rate, E , is defined bv 

n 

E » siny * V /g 
n 'a a 

where y is the aerodynamic flight-path angle and V is the true airspeed. One 
a a 

such schedule is illustrated in Figure 1, where circled points indicate data storage 
poirts. The reference controls obtained from these tables were derived subject 
to the maneuver margin requirements and engine RPM limits specified for STOLAND. 
The values of engirn RPM are also corrected to compensate for deviations of the 
ambient temperature from standard values. 

I'pon engagement of the 'Capture' mode (entered by pushing the Reference 
Flight Path button) the system synthesizes a minimum-fuel reference trajectory 
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from the current aircraft position lo touchdown. The reference flight path con- 
sists of two parts as shown in the example in f igure 2 . The fixed route is specified 
b\ input waypoints and the final segment is always a straight, constant -speed 
glide slope. Note that althou :h the trajectory is computed to touchdown as shown 
in the figure, the minimum-fuel guidance system will be disengaged during the final 
seg nent at the preset altitude sufficient to allow manual landing or go-around. 

The initial portion of the trajectory, shown in dashed lines, is the 'Capture' 
trajectory which transfers the aircraft from its current position to the capture 
waypoint specified by the pilot. The capture trajectory is resynthesized con- 
tinously as the aircraft moves until 'Track' mode is engaged (entered by pushing 
Fl'LLALTO or HORNAV button). At that time die reference trajectory is fixed 
and the aircraft is automatically controlled to the reference. There is no flight- 
director guided manual mode. 

The optimum guidance system can be initiated from any initial aircraft position 
heading or airspeed in the terminal area where a capture waypoint on a stored re- 
ference trajectory and a loiter speed is specified by the pilot. The system will 
attempt to synthesize a trajectory to the runway and provides information on the 
MFD display related to no-capture conditions. When capture is feasible, the pre- 
dicted time of arrival at the runway is displayed. 

In the 'Capture' mode, the pilot can use speed control or path streching to con- 
trol the approximate time of arrival. The trajectory is compensated for pre- 
dicted wind conditions entered by the pilot via the keyboard as a function of altitude. 
The display techniques and two dimensional capture trajectory generation method 
was taken fiom previous work (Hefcrence 3). 
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Capture Houte 
(resynthesized continually) 
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Pilot Operation 


The wind profile and runway temperature are entered through the keyboard 
a* follows 

RWT *xx xx • runway temperature, °C 

EW_-yy.xx xx* Wind magnitude, knots 

yy • Wind he ad in? degrees /10 


The third letter in E\\_ denote* altitude 
as follows: 


MNEMONIC ALTITUDE , FT 


EWA 

EWE 

EWC 

EWD 

EWE 

EWF 

EWG 


0 

250 

500 

1000 

2000 

4000 

8000 


These inputs will normally be made prior to the flight but can be changed 
during flight. 

When the aircraft reaches cruise altitude: 

AUTO FLAPS: ON 

AUTO ON 

AUTO THKOTTLE ON 

KEYBOARD : WGT « A/C WEIGHT, VBS. 

LAS * LANDING AlRSPfcLD. KNOTS 

( DEFAULT * STOLAND MINIMUM LAS) 

AIRSPEED HOLD i 

ALTITUDE HOLD J Optional variations in these parameters 

HEADING HOLD I will cause initial errors on track 

KEYBOARD LTS * LOITER SPEED, KNOTS 

(default value 140) 


-6 



PI SH TAC'AN Enable navigation system 

PISH KIM There is onh one reference flight path 

PISH REFP 


When Hi EP is puahed the system la in the ’Capture’ mode, and the fixed por- 
tion of the reference fltght path appear* on the MFD with the triangle which In- 
dicate » current aircraft poaition and heading. At th is time the capture waypoint, 
WP1 . muat be entered b\ keyboard. Th*- default value of WPT is zero and the 
guidance system will indicate no capture until a waypoint is chosen. WPT is set 
to zero whenever the REFP is turned on and must be keyboard entered. 

When a waypoint is selected the system predicts the aircraft poaition 15 
second* in thi future, assuming straight constant-speed flight, and attempts to 
construct a horizonttd flight path from that point to the capture waypoint. Fail- 
use to find a flight path results In a "NO CAP HOP." message on the MFD. This 
can be corrected bv changing aircraft position and heading or the capture way 
point. Any waypoint except the last one (touchdown) may be captured. A 15 
second prediction is made using the current aircr~.it speed, heading and flight 
path angle at the beginning of each capture trajectory to allow time for svnthesis 
computation, and variations in those parameters during the cycle will cause 
initial errors when tht track mode is engaged. 

When the horizontal capture path has been found, see Figure 3, a minimum- 
fuel speed-altitude profile is generated along the combined capture and fixed 
routines from aircraft to touchdown by integrating the point-mass equations of 
motion. Beginning at point P 0 the aircraft immediately accelerates or decelerates. 

m 

at constant altitude to the "loiter" speed, LTS, input by keyboard; airspeed is then 
held constant until the change to the capture waypoint speed must be initiated. 
Figure 3 provides a graphical representation of the speed-altitude profile to the 
capture waypoint. 
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Figure 3. Typical Speed - Altitude l*r**fi1e in Capture 


Each time the capture trajectory is resynthesized a new time of arrival, 
TOA,(at touchdown) is computed and displayed on the MFD. Hence, coarse 
adjustments to the YOA may be made by maneuvering the aircraft or changing 
the loiter speed or capture waypoint. When a satisfactory trajectory has been 
achieved the pilot engages the track mode by pushing either HORNAV or Ft’LL- 
Al"IO. The system may be returned to the predict mode any time by turning 
off F I'LL ALTO and HORNAV. Note that invoking one of the Sperry Autopilot 
modes while the track mode is engaged destroys the reference trajectory and 
will result in spurious autopilot response. Should another automatic mode be 
invoked Inadvertently during track, RE FI’ should be turned off and the procedure 
repeated. 

During track mode operation, when the aircraft reaches a preset altitude 
above the runway (presently 300 ft.) the minimum-fuel guidance systen. disconnects 
automatically leaving the STOLAND system in the standard control wheel steering 
(CWS) mode, and the autolhrottle and auto switches in the "on” position. The 
system can also lx- disengaged by moving the control column or pressing the 
autopilot disconnect button, in these cases the auto and the autothrottle switches 
are also turned off. 

The following keyboard entered parameters are related to system control 
along the reference trajectory. 

Trajectory synthesis parameters: 

limits 

GNA > - V /g ratio 0-100 

a 

GNB E margin for control 100-500 

Track mode operation parameters: 

GNC Throttle rate trip level 25-500 

GNI) RPM min bias 0-500 
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ONE 

Max deceleration in RFP 

limits 

0-200 

GNF 

Nozzle monitor trip level 

25-500 

GNG 

A alt feedback on pitch •lOOfps/V 

■ 

0-200 

GNU 

a 

A alt feedback on pitch *100fpgA’ 

* c 

0-200 

GM 

0 pitch rate gain 

25-500 

GNJ 

Nozzle feedback control factor 

0-200 


Control Logic 

Reference flaps are determined bv the optimization procedure at each 
point along the reference flight path. The flap command is the reference value 
limited to the placard value as determined from the current acU.M airspeed 
and also constrained to prevent retraction. Some combinations of capture way- 
point and loiter speed may call for extension of the reference flaps followed by 
retraction after reaching the capture waypoint. This will result in the actual 
flap setting larger than the reference value and cause slowly decaying transient 
er*ors in speed and altitude. The transient is repeated at each change in command 
until the reference flaps increase .again. 

Figure 4 shows the flight envelope divided according to the use of nozzle ar.d 

throttle for contrc 1. Over most of the flight envelope nozzle is used only when 

the throttle is at its lower limit and further reduction in E is called for. (For 

this purpose settings of 45 or more it is actually the reference value.) However, 

for large reference flap settings E is increased bv simulataneously increasing 

n 

the throttle and retracting the nozzle. 

The standard system stall prevention and RPM limiting routines are operative. 
However, the engine HPM’s commanded by the minimum-fuel system are generally 
within the Sperry limits and no noticable effects are seen from either RPM 
limiting or stall prevention duiing normal operation. The maximum RPM algorithm 
has been modified as required for revised engine data. 
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The standard airspeed limit no longer has a direct connection to 
the throttle servo. The "Minimum Speed" and "Haise Flaps" messages are 
still displayed when the standard limits are encountered. The "minimum speed" 
message is usually on In cases where flaps are extended and there is a large 
negative flight -path angle. In this case the standard minimum, includ'ng bant 
angle protection is 3 or 4 knots higher than the reference airspeed; however, 
the reference conditions always meet the maneuver margin requirements. 
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n AVIONICS SOFTWARE DESCRIPTION 


Executive for the Airborne System 

The program has been modularized to separate the capture trajectory 
synthesis segment of the system from the time constrained guidance and con- 
trol segment. The revised avionics software utilizes existing navigation, inner 
loop autopilot and many automatic modes with only miuor modifications. The 
basic capabilities removed for this application are the preflight test module and 
the reference flight path automatic guidance and control mode. 

To provide the required real time for the Optimum profile guidance and 
control algorithms a new executive was Incorporated in the basic program. The 
executive used has been tested previously in a flight environment in the Twin 
Otter Kalman Filter Airborne Program (Reference 4». The executive was de- 
signed to execute foreground and 4 levels of background. Foreground, a 10 Hz 
control module background and a far background synthesis module have been 
used for this application. To accommodate the levels of background, the execu- 
tive is designed to save the restore interrupt locations, register values and 
library function parameters for each level of background. The library function 
parameters allow the SINCOS, ARCTAN and SQROOT to be re-entrant at any 
level. A macro-flow chart of the executive is presented (Figure 5). At pro- 
gram start initialization logic is executed. The interrupts are then enabled and 
the far background logic is entered. 

The far background logic generates a continually updated capture trajectory 
to the specified capture waypoint and synthesizes this trajectory. During 'Track' 
mode, the far background module performs derivative calculations for the for- 
ward integration. The system performs these operations recursively using 
time not required to service foreground or other levels of background. The pro- 
gram performs far background- calculations until an interrupt occurs. The 
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MACHO FLOW CHART 


ADVANCED GUIDANCE AND CONTHOL EXECUTIVE 


I entry 


I 


[ INITIALIZATION 


Jl 


3 


fENABLE 1kHz INTERRUPT 1 
DTG4D flNffl A TE FAR BACKGROUND 1 


1kHz INTEHHUPT 


] 


TIME TO START FOREGROUND 
(20Hz) 


no 


-» i RETURN TO EXECUTION 


SAVE INTERRUPT RECOVERY 

IK PU SH rlK'WNl £1A£L 


£ 


I EXECUTE FOREGROUND 


l 


I 


TIME TO START 10Hz 

no 

RESTORE INTERRUPT 

BACKGROUND 


RECOVERY FROM STACK 


(ToHz OVERFLOW 

■■■ifiKra 


AGEXC | EXECUTE 1 OH 7 BACKGROUND 
»' ■ — 1 " 


1 yCS 

I 

STORE OVERFLOW 
LOCATION 

J * 



1 00 ft 

RESTORE INTERRUPT 

J * 

RECOVERY FROM STACK 

: .... _r 


STORE OVERFLOW 


LOCATION 


TIME TO START SECOND 
BACKGROUND 


| overflow" 


not used 
for this 
applica- 
tion 


I EXECUTE BACKGROUND ~| 


E 


TIME TO START THIRD 
BACKGROUND 


no 


RESTORE INTERRUPT 
RECOVERY FROM STACK 


("OVERFLOW 


yes 


I 


STORE OVERFLOW 



(location 

t 

DTG4I) 1 EXECUTE FAR 

BACKGROUND | 


Figure 5. Advanced Guidance and Control Executive 
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system requires several sources of interrupt, however, only the interrupts 
triggered by the 1kHz internal clock of the 1&19A computer are presented. A 
clock interrupt causes transfer to a location where a counter is decremented 
and tested to determine if it is time to initiate the 20 Hz logic (foreground) calcula- 
tions. If it is time to start the 20 Hz logic, then the necessary quantities for 
recovery from interrupt (registers, location at interrupt etc.) are saved in a 
software push down stack. The program then executes the foreground Augmentor 
W ing navigation guidance, control and display software. In foreground operation, 
interface outer-loop control and display parameters are transfered to the standard 
system from the Advanced Guidance and Control modules as required. Next, 
the logic determines if the 10 Hz calculations (every other cycle) should In- 
initiated. If not the executive branches to the saved location utilizing the push- 
down stack, Fach level of background is executed in a similar manner where at 
least some part of the far background logic is executed on each cycle. For this 
application, the control parameter calculations and the associated real time for- 
ward Integration of the trajectory’ are performed in the 10 Hz background segment. 


Optimum Profile Guidance and Control Module Structure 

A functional description of the software system is provided in this section. 
Functional diagrams of all software modules are presented for the synthesis, 
forward integration, control and display modules developed for this application. 

The figures provided a description of the general structure of the program including 
the subroutine calling sequence where the routines associated with the various 
modules are so indicated. Software modules in the standard STOLAND software that 
were modified to provide the required interface are specified in each structure diagram. 
Following the structure diagrams a subroutine glossary is presented providing a 
functional description of the purpose of each subroutine. 
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OPTIMUM PROFILE GUIDANCE CA PTURE/DERIVATIVE MODULE 
(EXECUTED IN FAR BACKGROUND) 
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* sub-programs provided In basic avionics software 
( ) Interfacing sub-programs 



OPTIMUM PROFILE GUIDANCE MFD DISPLAY MODULE 
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* *ub -programs provided In butc avionic* scflware 
( ) intrrf.irinp sub-propramn 


Subroutine Functional Inscription* 


AGEXC 
ARCKI'T 
AHCS1N 
AT ANA 
ATM FLT 
AT MINT 


BLLM 

CAPTHJ 

CNTHL 

CNTHLI 

CTURN 

DSCADV 

DTG4D 


Ewnnis 


Optimum profile guidance lot) MS background executiv. Driver for 
the forward integration and control sub-program* 

Loads and convert* aircraft parameters for the autopilot control 
sub-program. 

Approximate arc sine routine. Input sine at BO, output angle in 
radians. 

Approximate arc tangent. Input tangent at BO, output angle in 
radians. 

Computes atmospheric parameter* as a function of altitude and 
temperature for reference trajectory calculations. 

Computes Initialization parameters used in atmospheric calcula- 
tions including temperature correction coefficient for landing 
airspeed calculation. 

Function to limit between input maximum and minimum values. 

Generates the capture trajectory for the MFD display. 

Calculates control parameters for the outerloop autopilot using 
the advanced guidance synthesis ard current aircraft configuration. 

Provides first entry to CNTRL at Track mode initialization. 

Computes turn points for capture trajectory MFD display. 

Descale angle for display function. 

Advanced guidance far background executive. Driver for the 
trajectory synthesis during capture mode and derivative calcula- 
tions for forward integration in Track mode. 

Provides the capability to display the values of EWA - EWG. 
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EW DLOD 

FLMX 

FLMXC 

FLNN 

FODDIS 

FODLOD 

INTEG 

INTP.PL 

JTHRTL 

Kimicx 

LTSD1S 
LTSLOP 
MI NSEC 
MSC4D 
NEWPSI 

NOMTRJ 


Provide* the capalility to load ar estimated wind profile aa a 
function of altitud* via keyboard entry (EWA - EVtC) 

function uacd to calculate reference value of flap* aa a function 
of air speed. 

Performj same function as FLMX for the 100 MS module. 

'imputes a throttle fuel rate factor aa a function of throttle and 
tempe rature. 

Displays keyboard FOD mode option. FOD not currently activated. 

Provides the capability to load the FOD mode option via key- 
board. FOD not currently activated. 

Performs the altitude - speed profile integration used in the 
trajectory synthesis. 

Performs a linear interpolation for a function of a single variable. 

Provides the interface between the advanced guidance control 
module and the inner loop autopilot for throttle and nozzle. 

Performs X* /3 trij. .lomelric expansion coefficient calculation. 

Displays capture loiter speed via keyboard. 

Loads capture loiter speed (LTS). 

Formats time parameters for MFD display (TWPT, TER). 

Optimum profile guidance MFD display executive. 

Computes the turn parameters for the 2D capture trajectory. Used 
In initialization to compute turn parameters for the fixed trajectory. 

Performs derivative calculations for the forward integration in 
Track mode. 
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NOM I HJ 

NT MOD 
NUKBHD 

PAULK 

PIMOD 

PHEFFP 

KWTDIF 
HWTLOD 
SC HCH1 
SEHCH1 
SINCS1 

SI NX 

SQItTA 
SI PINT 

SI ItTEH 
STKWYD 
TANX 
T DA YET 


Performs the real lime forward integration of (he reference 
trajccloi) in track inode. 

Performs angle modulation for NOMTK2. 

Provides keyboard enir)' of two parameters for wind profile 
entry (heading and speed) 

Pitch arming routine. 

Performs angle modulation for far background module. 

Performs pilch engage for the advanced guidance reference 
flight path mode. 

Displays runway temperature via keyboard. 

Loads runway temperature via keyboard. 

Performs same operation at* SLHCH1 for the 100 MS module. 

Computes algebraic function used in linear Interpolation. 

Interface to re-entrant SINCOS routine with index register 
protection. 

Approximate sine routine. Input angle in radians at BO, output 
at BO. 

Interface to re-e trant SQHT routine with index register protection. 

Computes derivatives and integral « the reference trajectory 
a single step. 

Competes time error MFD display. 

Displays capture waypoint parameters on MFD. 

Approximate tangent routine. 

Computes time of day for predicted arrival at runway. Used for 
MTD display. 
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THAN I M 
TKKSE ! 
TST 
TWOD 
VHTSYN 

W1NMOD 

WNDCAL 

WPBGND 

WPENTR 

WPSPD 

WYPTDS 


Performs TACAN lo MODU.S transition. 

Stores capture waypoint parameters into MM) waypoint array. 
Computes the 2D capture trajectory. 

Computes the 2D fixed trajectory. 

Controls the compulation of the altitude - speed profile for the 
complete reference trajectory and loads the command table. 

Extracts wind parameters from estimated wind table for reference 

trajectory. 

Generates estimated wind table from keyboard entered data. 

Transfers reference flight path parameters to advanced guidance 
internal tables and performs Initial calculations for MFD display. 

Transfers waypoint coordinates. 

Converts waypoint speeds to internal units and applies landing 
airspeed constraints where required. 

Displays lime of arrival and time error on MFD. 
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for* 1 Sumnun of th« Airborne System 


rhc* table provide^ below presents a core aumman of the complete system 
including the- revised standard software and the Optimum profile guidance and 
control modules. The total memory storage used for the flight software is 
28082 words. Total system capability is 32K words of which 2000 words an re- 
served for utility purposes. 


The first seven sections represent the revised standard system and the 
fin.d section! represent the new capability implemented. 


Kevised Standard System 


Section 

Locations 

Si/e 

1 

30000 

. 

3377(1 

2047 

2 

34000 

- 

37750 

2023 

3 

200 

- 

1241 

546 

4 

10000 

- 

162(> 1 

3314 

5 

20000 

- 

27300 

3777 

6 

40000 

- 

47743 

40(1 * 

i 

60000 

- 

67102 

36. ’>1 


Advanced 

Guidance 

ani- 

Section 

Locations 

Size 

(4 

124(1 

• 

2747 

834 

9 

16400 

- 

17630 

665 

10 

27350 

- 

27714 

9*>q 

11 

50000 

- 

57552 

3947 

12 

72000 

- 

77642 

2979 


Description 

P'tnction lib ran - Data Art a 

Protected Data 

Executive 

Autopilot 

Data Interface 

Kaviga.ion - Displays 

MFD Display 

rol Modules 

Description 

Capture Modules 
Autopilot - Control 
Keyboard Entry 

Data - MFD Interface - Synthesis 
Derivative Integration 
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Ill FORTRAN SIMULATION PROGRAM (AIG4D) 


Structure and Design 

The final section of this report represents a current update ami extension to 
AMA Report No. 7U-iy performed under contract NASli-yjlO (Reference 3). 

Modular Str u cture The purpose of this Fortran simulation program is to 
assess the performance of a fuel efficient onboard guidance system for powered 
lift SI OL aircraft. The simulation program, resident on the Ames TSS:i60/G7 
facility, provides a test-bed for the Optimum Profile Guidance and Control Air- 
borne System. The Optimum Profile Guidance Simulation Program, Al’G4D, is 
a combination of a fast -time 3D guidance simulation and a 4D trajectory synthesis 
program developed by NASA/ARC. The program has been modularized to sepa- 
rate the onboard portion of the simulation as presented in Figure y. The sub- 
routine calling sequence is presented in Figure 10. This figure describes the 
general structure of the program where the routines associated with the various 
modules are so indicated. The display and real time data entry software provided 
in the avionics software has not been simulated in the FORTRAN Program. 

Following this discussion a subroutine glossary is presented giving a brief 
description of the purpose of each subroutine. Finally, sub-program and common 
cross-reference tables for the complete program are presented in this section 
In this implementation, the common arrays have been grouped so that each labelled 
common is associated with a particular module. The labelled common arrays 
required by the flight module are indicated in the Common Array Size Require- 
ment Table. 


AIG4D MODULAR ST RUCTURE 



Figure 9. 


AUG4D Modular Structure 










SUBROUTINE CROSS REFERENCE TABLE 


nanc subprograms referencing renbcr 

ACRSIN RAIN 


ASEXC ACRSIR 

ARC 1 M T ACRSIR 

ARCRFT CNTRL 

ARBC1 ATHFLT ATNSPM RATEUF 

ATNINT 1ST 

ATNSPM ARCRFT 

ATRSTP CNTRL INTER NONTRJ TST 

BUR CNTRL 

CHKSTP RAIN 

CNTRL AGEXC 

CNTRL I AGEXC 

E MGDAT ARCRFT 

FL IN ENGDAT RATE WF 

FLNX ARCRFT CNTRL STPINT 

FLNN nortrj 

F 3D ARCRFT ENGDAT 

I N I T RAIN 

1NTE6 VMTSYN 

NE UPS 1 TST TUOD 

NONTRJ AGEXC 

NONTR2 AGEXC 

PIP10D ACRSIR ARCRFT CNTRL NEUPSI TST TUOO 
RATEUF INTEG 

REFUF NONTRJ 

SERCH1 ARCRFT ENGDAT FLNN RATEUF STPINT TST 
UNDFLT 
SGN NEUPSI 

STPINT ARCRFT CHKSTP INTE6 NONTRJ 

TREAD INI T 

TST NAIN 

TUOD INIT 

VHTSTN VSYN 

VSTN RAIN 

UINDAT INIT 

UINNOD INTE6 NONTRJ 

UNDFLT ARCRFT 


ENTRY POINT SUNNARY 


ENTRY 

SUB-PR06RAR 

ARCINT 

ARCRFT 

ATNINT 

ATRSTP 

CNTRL 1 

CNTRL 

NONTR? 

NONTRJ 


Figure 11. Subroutine Cross Reference Table 


UINNOD 
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CONAON CROSS REFERENCE TABLE 


NAME 

IUIPROIRANS REFEtENC INS NEARER 



ACIATA 

ARCRFT 

BLR BT 

ENIDAT 

IAIT 




AC F L T 

AC RSI 4 

ARCRF T 

crtrl 

INI T 

NONTRJ 

IMBFLT 


ACREF 

ACRSIR 

AOEXC 

ARCRFT 

CNTRL 

IAIT 

NOATRJ 


II 

INI T 

INTES 

NOATRJ 

TST 

TUOD 

VSTN 

U1NA0D 

1 1 A 

ARCRFT 

1N1T 

VSTN 

undflt 




12 

ARCRF I 

BLK dt 

CHKSTR 

CNTRL 

INI T 

INTEI 

HAIN 


NOATRJ 

STRUT 






13 

ARCRFT 

INI T 

INTEI 

MAIN 

NOATRJ 

VHTSTN 


14 

ARCRFT 

chrstp 

INI T 

INTE6 

NOATRJ 

STPINT 


CAFIT 

ATNFlT 

BLK DT 

CHRSTP 

FLAN 

INTEC 

NONTRJ 

RATEUF 


STRINT 

TST 

TUOD 





CAFLTA 

ATNSTR 

CNTRL 

REFUF 





CONTRL 

ACRSIN 

ARCRFT 

BLK DT 

CNTRL 

HAIR 



01 

INI T 

MAIN 

TST 

TUOD 

VHTSTN 

VSTN 


ENDATA 

ARCRFT 

ATNSRH 

BLK DT 

CNTRL 

EN6DAT 



INOUT 

ACRSIN 

AGE XC 

ARCRFT 

ataflt 

ATASPH 

ATASTP 

BLK DT 


CMKSTP 

CNTRi 

E4GDAT 

INI T 

INTEG 

HA IN 

NOATRJ 


STMAT 

TST 

”<UOD 

VHTSTN 




i Aici 

ACRSIN 

CNTRL 

NOATRJ 





IATOI 

INTEL- 

VHTSTN 






STOL 

CHK5TR 

INI T 

HAIN 

STRINT 




STM 

ACRSIN 

ARCRFT 

CHRSTP 

INTE6 

NONTRJ 

STRINT 


STA 

ARCRFT 

VSTA 

INI T 

MAIN 

NOATRJ 

TST 

TUOD 

VHTSTN 


Fipure 12. Common Cross Reference Table 
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Subroutine Glossary 


MAIN 

AC RSI M 

AGEXC 

ARCRFT 

ATMFLT 

ATMSPH 

ATMST P 

BUM 

CHKSTP 

CNTRL 

ENG DAI 

FUN 

FLNN 

FLMX 

F3D 

IMT 

INTEG 


NEWPS1 


NOMTRJ 
PI MOD 


Provides the overall control for the execution of a case 
Executive for the aircraft simulation 

Executive for the forward integration and control sub-programs 
Aircraft simulator 

Evaluates the atmospheric model for the flight module 

Evaluates the atmospheric model for the simulation module 

Evaluates the atmospheric model for the onboard module 

Function to limit between input minimum and maximum values 

Software verification module to evaluate derivative calculation? 

Control law routine, generates commands to fly the aircraft 
along the reference 

Engine model routine used by the aircraft simulation module 

Algebraic function used for linear interpolation 

Computes fuel rale as a function of throttle and temperature 

Function used to calculate flaps as a function of speed 

Algebraic function used for linear interpolation of a function 
of two or three variables 

Controls the input of data for the preflight data load 

Determines the altitude-speed profile used in the flight module 
trajectory synthesis 

Computes the turn parameters for the 2D horizontal trajectory 

Reproduces the synthesised reference trajectory in short 
segments using the command table 

Function to reduce an angle to its principal value t n radians 
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KATKW F 
RFFW F 
SERCH1 
SGN 
ST PINT 

TREAD 

TST 

TWOD 

VHTSYN 

VSYN 

WIND AT 

WNDFLT 

WINMOD 


Compute!! fuel rate for thi aircraft moduli 

Compute!* fuel rate for this onboard module 

Computes an algebraic function used in linear interpolation 

Function to return the sign of a real variable 

Computes the derivatives and integrates the aircraft trajectory 
a single step 

Controls program input of the list, drag and elevator used in 
the aircraft simulation module 

Computes the 2D capture trajectory 

Computes the 2D fixed trajectory for the pre-flight data load 
module 

Controls the computation of the altitude-speed profile for the 
complete trajectory and loads the command table 

Controls the compulation of the minimum, maximum and 
nominal time trajectory’s and loads the 4D 'ables 

Controls the input of the wind tables 

Performs wind calculations for the aircraft simulation module 
Performs wind calculations for the flight module 
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Program Input 


In this faction, the Inputs, and program options available through appro* 
priate input selection are described. The inputs ari conveniently divided into 
two groups or blocks, the pre-flight data load block read on logical units 8 and 
the standard input block read on unit 5. The pre-flight input block is composed 
principally of tabular data Input via fixed formats. The pre-flight data block 
is not varied from case to case. The standard input block is read using the 
namelist feature of the Fortran language. The namelist is AUG4. 

P re- Flight Data Load Input 
Data Description 

HEADER (20) Any alpha-numeric information 

FORMAT (2014) 

Wind Tables 

ALTW (12) Values of altitude for the INWIND table 


FORMAT 2(GF10. 2) 

INWIND (48) 1-12 Estimated wind speeds knots 

13-24 Estimated w ind headings deg 

25-3G Actual wind speeds knots 

37-48 Actual wind headings deg 

FORMAT 4(1215) 


The "actual" wind parameters are used in the 
simulation module. 

Waypoint Table 

NWPI Number of input waypoints 

FOR MAT(14) 
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XWP 

Vlaypoint x coordinate 

ft 

YWP 

Waypoint > coordinate 

ft 

ZWP 

W aypotnt z coordinate 

ft 

R 

Waypoint turning radius 

ft 

VNOM 

Nominal waypoint speed 

knots 

V MAX 

Maximum waypoint speed 

knots 

V MIN 

Minimum waypoint speed 

knots 

FORMAT N\YP1(10X,3F10. 1.4F10.4) 
Energy-rate Schedule '1 able 

M1NMAX Throttle variation index limit G 

LVMAX 

Speed variation index limit 

G 

LHMAX 

Altitude variation index limit 

!. JJ 

LVYMAX 

Weight variation index limit 

* 2 

MAXCOF 

Coefficient matrix index 

s 6 


i'he matrix index is specified from the 
following table: 

1. Speed 

2. Nozzle 

3. Normalized energy rate, E^ 

4. Interpolation coefficients 

5. Angle of attack 

G. lnterploation coeffients 

FORMAT (514) 

HSET(LHMAX) Altitude table for energy-rate schedule ft 

table 

FORMAT (2E14.7) 
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WTSET(LWMAX) Weight table for force schedule table |l>t» 

FORMA! (2E14.7) 

VU'HA (M1NMAX, LVMAX, LHMAX, LWMAX, MAXCOF) 

Force schedule matrices used to obtain 
nozzle, and angle of attack. 

FORMAT (3E14.7) 


Lift and Drag Tables 

NF Flaps variation index limit 

NCJ Cold thrust index limit 

NALFNV Angle of attack index limit 

FORMAT (3I5> 

CLP (NF, NCJ, NALFNV) 

Lift coefficient table 

CDP (KF.NCH, NALFNV) 

Drag coefficient table 

DEP (NF, NCJ, NALFNV) 

Elevator table deg 

Namelist AUG 4 


DACTS 

Time step size used in the integration 
of the aircraft simulation equations 

sec 

. 1 

DEDS 

Distance step size for the integration 
of the synthesized trajectory 

ft 

100. 
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DEDTK 

Gain on A> used in throttle control 
equation 


.17 

DEINTK 

Gain on integral of error in speed used 
in the throttle control equation 


. 13E-2 

DKX 

Gain on error in altitude in control 
law for throttle control equation 


.035 

DLT 

Time offset for capture initialization 

sec 

10. 

ELM AX 

Maximum limit on elevator 

deg 

♦15. 

El.MIN 

Minimum limit on elevator 

deg 

•23. 

ENAK 

Gain on altitude error in nozzle control 


.5 

ENGK 

Gain on flight path error in nozzle control 


3. 

ENSK 

Gain in integral speed error in nozzle 
control 


.5 

ENUK 

Gain on nozzle in aircraft simulation 
model - servo gain 


1. 

ENt'LK 

Nozzle deployment lead time ratio 


.5 

ENUMAX 

Maximum limit on nozzle 

deg 

104. 

ENl'MIN 

Minimum limit on nozzle 

deg 

6. 

ENURL 

Rate limit on nozzle 

deg sec 

5. 

ENVK 

Gain on speed error in nozzle control 


5. 

EPSLN 

Control factor for E 


1. 

ETMAX 

Maximum limit on E 

n 


1. 

ETMIN 

Minimum limit on E 

n 


-1. 

FAK 

Gain on flaps in aircraft simulation 
model - servo gain 

deg 

.5 
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FDMAX 

Matt limit on flaps 

deg 

sec 3.5 

FLK 

Flap deployment lead time ratio 


.5 

FLMX2 

Flap limit when throttle us less than 

8». 5', 

deg 

43. 5 

GAMK 

Gain on altitude in flight path angle 
control 


.1 

GLK 

Lead time ratio for flight path angle 
changes 


.5 

HAC 

Initial aircraft heading 

rad 

* 

I4D 

4D control flag. Set to zero only 3D 
synthesis is performed 


0 

1C)PT 

Simulation control flag. Set to one 
generates reference and simulated actual 
using ACRSIM. If set to two only the 
reference trajeclor) is synthesized 


1 

1SCALL 

Provides option to generate energy-rate 
schi le table file for use in onboard 
program 


0 

IPRNT 

Provides detailed print of derivative 
parameters 


2 

ISTND 

Provides stand alone derivative cal- 
culation capability 


0 

KKOUNT 

Print step bypass counter in ACRSIM 


20 

KMODE 

Aircraft simulation flag. Set to one 
the aircraft is simulated using ARCRFT. 
Set to two the simulation is bypassed 


1 

KPR1NT 

Trajectory synthesis print control. 

Cc! to zero generates most complete output 


0 

KSTEP 

Number of integration setps using constant 
derivatives in generating the reference 
trajectory 


10 

NCAP 

Capture waypoint index 


1 
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PLIM 

Roll rati limit 

deg set 

10. 

PLK 

Lead time ratio for bank angle change* 


• 7 

PHPSK 

Gain on cross track velocity error in 
roll equation 


1 . 

PH VK 

Gain on cioss track position error in 
roll equation 


.002 

PK 

Gain on roll in aircraft simulation model 


1 . 

QLIM 

Rate limit for pitch 

deg ^tec 

5. 

HQ 

Fraction of energy rate allowed for 
control purposes 


.1 

RWT 

Runway temperature 


15. 

5GMKN 

Minimum limit on sine of flight path 
angle 


-. 130- 1 

SGMXX 

Maximum limit on sine of flight path 
angle 


. 130t*l 

TAN FI 

Tangent of maximum bank angle 


.46031 

TBLSR 

Energy rate schedule table temperature 
lapse rate 

deg /ft 

0019*12 

THACK 

Gain on throttle in aircraft simulation 
servo model 


1 . 

THARL 

Throttle rate limit (not activated) 

/sec 

o 
•» • 

THDK 

Gain on altitude error feedback in pitch 
control 


10. 

THEK 

Gain in energy error in throttle equation 
in control law 


.023 

THOGK 

Not activated 


0 . 
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TIHVK 

Rate limit time constant for pitch in 
aircraft simulation model 


10. 

T1MEDA 

l)ela> time for trajectory synthesis 
(not activated) 


0 . 

TTBSL 

Standard sea level temperature 

°K 

2ms. 3 

TMPCR 

Cruise temperature 

°C 

15 

VAKC 

Initial aircraft speed 

knots 

• 

\ GMIN 

Minimum limit on true uirspeed rate 

Gravity normalized 

•ec 

*.03 

VGMX 

Maximum limit on true airspeed rate. 
Gravity normalized 

tec 

.05 

VLN 

Nominal aircrai' loitering speed 

knots 

133. 

VMAXN 

Maximum aircraft loitering speed 

knots 

160. 

VMINN 

Minimum aircraft loitering speed 

knots 

130. 

WGT 

Aircraft weight 

lbs 

3H000. 

XAC 

Initial aircraft x coordinate 

ft 

4 

YAC 

Initial aircraft y coordinate 

ft 

4 

ZAC 

Initial aircraft z coordinate 

ft 

4 

PHLIM 

Maximum bank angle allowed in 
control 

deg 

30. 

DH1NTK 

Gain in integral altitude feedback 
in throttle control 


.001 . 

THIK 

Gain in integral altitude feedback 
in pitch control 


. 0033 

EN1IIK 

Gain in integral altitude feedback 
in nozzle control 


.04 

gnit; 

Nozzle control feedback factor 


1 . 


• no default is provided for these parameters. They must be present in the 
NAMELIST input. 
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Kamel if] si PM 


Required Input parameter* arc indicated by an •. 

ALT Aircraft altitude ft 

CEPS Altitude. 'speed change ratio 

COSPH1 Cosine of bank angle 

E DT MAX Maximum change in energy rate. 

Gravity normalized 

EDTM1N Minimum change in energy rate. 

Gravity normalized 

EPS Control factor for energy rate 

RTCFE Temperature correction factor 

RR Fraction of energy rale reserved for 

control pu eposes 

SGMX Sine of maximum aerodynamic flight 

path angle 

SGMIK Sine of minimum aerodynamic flight 

path angle 

VLAS Indicated airspeed ft /sec 

YDGMX Maximum deceleration. 

Gravity normalized 

VDGMK Minimum deceleration. 

Gravity normalized 

VT True airspeed ft/sec 

WT Load factor lbs. 

AKW Partial of wind speed with respect to 

altitude 

KN1 RJ STPINT angle of attack computation flag 

KSTOL STOL mode indication 

ALPHA Angle of attack deg 
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ENT AC 

FAC 

SING AM 

THAC 

THM1N 

VADTG 

VIAST 

DEE 

FLMXJ 

FT MAX 

ETMIN 

IFRNT 


Nozzle setting deg 

Flap setting deg 

Sine of flight path angle 

Throttle setting r * 

Minimum throttle setting '» 

Nominal deceleration. 

Gravity normalized 

Speed criteria for flap deployment 

Total energy change on current waypoint 

Flap limit when throttle is less than deg 

H9. r> 

Maximum limit on energy rate. 

Gravity normalized 

Minimum limit on energy rate. 

G ravit\ normalized 

Print flag for additional output in ST PI NT 


♦ 






Deck Setup and Machine Requirements 

The AUG4D trajectory program and associated input data sets are 
stored at the NASA Ames Research Center TSS computing facility. The 
Fortran source decks are stored using standard line formal. The object 
modules are stored as a single job library presently named TEX. Following 
is a typical deck setup for executing the Al'GID using the loader facility' of 
the TSS operating system. 
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Program Execution 


AMES USYSLIB 

obtain standard atmosphere using system moduli 1 
ARDC1. 

JBLB TEX define job library for the AUG4D object modules 

DDEF FT08F001 , , AUGYDAT 

define inp’‘t unit for extremized force schedule table 

DDE F FTO'.FOOl , , AUG FIX 

define input unit for namelist AUG-4 input 

DDEF FT09F001, , FP4MDAT 

define input unit for pre-flight data load input 

DDEF FT10F001 , , STPNL 

define input unit for namelist STPNL 

DDEF FT11F001, .DTA1D 

define output unit for scaled force-schedule table used in 
onboard data load 

All program printed output is performed on unit l> and is printed on the 
installation defined standard output device. 

LOAD BLCKAW$:j> construct executable module 

LOAD MAINAWSS 

CALL MATNAW?* execute AUG4D 


The following tables present the total machine requirements of the Al t'. ID 
program with the current TSS implementation. Since the flight module is in- 
corporated in the Augmentor Wing STOLAND system, those associated sub- 
programs and commons are so indicated. 
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Suhro utinr S>/c Require nienlh 


MA1NAWX? 

BYTES 

3G1G 

FLT 

BLCKAWX? 

0 


ACRSIMXS 

2393G 


AGEXCX $ 

480 

* 

ARCRFTX? 

5888 


ATMFLTRS 

393 


ATMSPHR$ 

5.18 


ATMSTPXS 

500 

• 

BLIMXS 

280 

♦ 

CHKSTPX$ 

91G 


CNTRLXS 

2772 

♦ 

ENGDATRS 

1280 


FLINRS 

240 

* 

FLMXRS 

448 

* 

FLNNXS 

352 

♦ 

F3DRS 

8G4 


INTI'X? 

11784 


INTEGXS 

9520 

* 

NEWPSIR? 

1988 

+ 

NOMTRJX* 

8708 

* 

P1MODRS 

392 

* 

RATEWFRS 

588 

* 
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REKW X8 

512 

# 

SE HOURS 

510 

• 

SGNRS 

284 

• 

ST P1NTXS 

7476 

# 

TREAD3RS 

201(1 


TSTXS 

3692 

• 

TWODXS 

131(1 


VHTSYNXS 

4768 

» 

VSYNXS 

2492 

• 

WINDATRS 

237(1 


WINMODRS 

1284 

♦ 

WINFLTRS 

1160 



•Sub-programs required in the flight module. 


Common Arra> 

Size Requirements 


FLT 

LENGTH 

COMMON 

MODULE 

(WORDS) 

ACDATA 


4233 

ACFLT 


25 

ACREF 

* 

25 

B1 

4 

61 


( = 11213 words) 



B1A 

61 

H2 • 

ia 

b:< • 

984 

B4 • 

500 

CMKLT * 

41 

CMFLTA * 

3 

CONTRL * 

19 

D1 • 

192 

END AT A 

434 

INOl’T 

120 

INTCL * 

o 

INTG 1 • 

is 

STOL * 

4 

ST PI * 

27 

SYN ♦ 

178 


* 

Total data arrays required in the flight module. <2187 words) 



